Mott transition and integrable lattice models in two dimensions 
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We describe the two-dimensional Mott transition in a Hubbard-like model with nearest neighbors 
interactions based on a recent solution to the Zamolodchikov tetrahedron equation, which extends 
the notion of integrability to two-dimensional lattice systems. At the Mott transition, we find that 
the system is in a d-density wave or staggered flux phase that can be described by a double Chern 
Simons effective theory with symmetry su(2) 1 <gisu(2) 1 . The Mott transition is of topological nature, 
characterized by the emergence of vortices in antiferromagnetic arrays interacting strongly with the 
electric charges and an electric-magnetic duality. We also consider the effect of small doping on this 
theory and show that it leads to a quantum gas-liquid coexistence phase, which belongs to the Ising 
universality class and which is consistent with several experimental observations. 
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In spite of substantial advances in our theoretical un- 
derstanding of strongly correlated electron systems, sev- 
eral problems still continue to provide estimulating chal- 
lenges. One of the most interesting among these is the 
Mott transition, or metal-insulator transition driven by 
correlations. As early as in 1939, Mott argued that if the 
electron density in a metallic system was lowered enough, 
the Coulomb repulsion would dominate over the kinetic 
energy so that the system would undergo a transition to 
an insulating regime [l[ . From the experimental point of 
view, there exist several systems which display a Mott- 
type transition, such as vanadium oxide V2O3, several or- 
ganic conductors, some doped semiconductors, and even 
underdoped high Tc superconductors. Moreover, coexis- 
tence between phases of different densities has been ob- 
served in several experiments Q. From the theoretical 
point of view, finding solutions to even the simplest mod- 
els (such as the Hubbard model) is difficult, given the 
failure of pcrturbative approaches due to the narrow dif- 
ferences separating the localized regime of the electrons 
in the insulating phase and the itinerant one in the con- 
ducting state. 

As of today, there exist two basic approaches for study- 
ing this transition: one is the dynamical mean field the- 
ory (DMFT) met hod 0], valid in the limit of infinite di- 
mensions (or infinite coordination number), which maps 
the Hubbard model onto the impurity Anderson model, 
with the addition of a self-consistency condition. This 
framework neglects spatial correlations while retaining 
the on-site quantum ones. The second approach consists 
in finding analytic expressions for physical observables in 
integrable models exhibiting the Mott behavior, for ex- 
ample, by using the Bethe Ansatz or the bosonization 
methods [3] [H [a]. However, the main restriction of these 
models is that they are formulated in one spatial dimen- 
sion, unlike most system of experimental interest. The 



goal of the present article is to extend this approach to a 
two-dimensional (integrable) lattice model that exhibits 
the Mott transition and writ down an Efective fiel Theory 
to futher analyze the behavior af the model. 

Let us start by consider a system of spinless fermions 
on a square (two-dimensional) lattice with hamiltonian: 



H = --J2[TpHi + e u )e iA >"ilj(i) + h.c.} 
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-U^2p(i)p(i + e M ) 



(1) 



where i labels the lattice sites and e M are unit lattice 
vectors, t is the hopping parameter, U is the (constant) 
Coulomb potential, p(i) is the normal ordered charge 
density with respect to the half-filling ground state, 
p[i) =: ipjipi : —1/2 and is the abelian statistical 
gauge field, which after imposing Gauss' law constraint 
reads 



4.(0 = Et e ( fc ' - e ( fc ' 1 + ^Mv^ , (2) 

where Q(k,i) is the angle between the chosen direction 
i and an arbitrary one k on the lattice. Note that, for 
the one-dimensional case, the gauge field is irrelevant, in 
agreement with the fact that quantum statistics in one 
spatial dimension is arbitrary 7] and does not involve any 
physical gauge field. Using the two-dimensional Jordan- 
Wigner transformationQQ : 



s z = 



(3) 



where U 2 d(j) = cxp [i J2k^j ; j)tp k ipk],the hamilto- 
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nian ([ljbecomes that of a XX Z Heisenberg model 

Hxxz = £ S ! + S i S ]) + AS ! S ^ > ( 4 ) 

where we have rescaled the terms such that A = U/t. 
Following we define an interaction star as the set of 
points where the spins entering in an elementary interac- 
tion are localized, i.e., the central site and their nearest- 
neighbors in the XXZ model. The n-th interaction star 
has an energy Exxz([&]n) which depends on the spin 
configuration in the star and on the local Boltzmann 
weights W([cr] n ) . Therefore, the partition function takes 
the form Z = ^2 a Yl n W([a] n ), where the sum is taken 
over all possible configurations of the entire lattice. 

In two dimensional quantum systems and three dimen- 
sional statistical models the integrability is guaranteed by 
the existence of a set of mutual commuting layer-to-layer 
transfer matrices T mn (l,fi), which is tantamount to the 
existence of solutions of the so-called Zamolodchikov's 
tetrahedron equation(TE) [H : 



RabcRadeRbdf Reef — Reef Rbdf RadcRabc 



(5) 



where the operators Rijk define the mapping Rijk ■ 
Vi <8> Vj ® Vfe — ► Vi ® Vj ® Vfej and V n is the spin one- half 
representation space , such that their matrix elements are 
the Boltzmann weights of the vertex Rijk = 
the indices i,j,k label the interaction star). These can 
be rewritten as the LLLR — RLLL operator conditions, 
which express the associativity of the Zamolodchikov al- 
gebra: 



Ll2,aL\z,bL2Z,cRabc — -Ro6c-^23c-^13 ,c-^13 .1 



(6) 



where, for example, the operator Li 2a acts on V\®V2®F a , 
Vi, V2 are the auxiliary spaces and F a is the quantum 
space. If F a is the representation space of some alge- 
bra A, it is possible to interpret the operators Lij a as 
operator- valued matrices acting on V\ ® V2, and depend- 
ing 'parametrically' on the generators of the algebra A 
denoted by v a and, possibility, on some c-numbers de- 
noted by s a - L\2a = Li2{v a , s a ) . In this case, the equa- 
tion © can be expressed as a 'local Yang Baxter' equa- 
tion: 

£l2(v a ,S a )Li 3 (v 6 ,S & )£23(v c ,S c ) = 

L 23 {'v l c,s c )L 13 (V b ,s b )L la (V a ,s a ) . (7) 

The tetrahedron equation §5§ is highly non-trivial to 
solve, but recently a new solution to it has been found in 
111] , The solution is associated to the finite-dimensional 
highest-weight representations of the quantum affine al- 
gebras U q (sl(2)), displaying the three-dimensional struc- 
ture of quantum groups. It may be understood as a 
quantization of the spatial fluctuations of geometrical ex- 
tended objects, and we shall see that in our case that 



these may be reinterpreted as (discrete) charge density 
waves. For completeness, we now briefly review the new 
solution (for details see [ill). The solution is inspired in 
the geometry of transformations applied to an hexahe- 
dron (see fig. |T|): there are three independent angles on 
each face and nine angles to fix the spatial orientation of 
the hexahedron. Therefore, nine independent angles are 
needed to specify it. Let us consider the mapping, 
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23 



[c4>M> ( g ) 
where a, (3 and 7 are the angles of the j-th face, and 
the primed variables refer to the opposite faces. For each 




FIG. 1: Diagram of a hexahedron in which the angles (a, (5, 7) 
and side lenghts (l p ,l' p ,l q ,l' q ) for one face are indicated. 

quadrilateral face, say the 1', the relationship between 
opposite sides is given by (I' , I' Y = X(Ai)(l p , Iq) 1 , where 
X(„4i) is a matrix acting non trivialy on the face V , that 
depends on the angles on that face (Ai = (cti,/3ij7i))j 
and which for a circular lattice reads: 



X{Ax) 



ki a\ 
—ai ki 
1 



where fci = cosaisin/3i, a = cosaisin(ai + Pi) and 
a* = cosai sin(ai — (3\). In the general case, considering 
three faces , we have: 

(l' p , l' q J' r f = X pq (A 1 )X 9r (A> i )X r .(Aa)(l p , l g , l r f ■ (9) 

The same result (|9]) is obtained by using the op- 
posite faces with angles {A'±,A' 2 ,A'z). Therefore 
, it is easy to see that there exists a func- 
tional mapping given by X pq (Ai)X qr (A2)X rs (A 3 ) = 
Ri23(X rs (A 3 )X qr (A2)X pq (A 3 )). This relation could be 
considered as a 'gauge symmmetry'. It can be shown 
that the mapping R123 satisfies the functional tetrahe- 
dron equation (FTE), a 'classical version' of the tetra- 
hedron equation: 



k' 2 a'i = k 3 a\ — k\a 2 a 3 
a 2 — a 2 a 3 + k±k 3 a 2 
k' 2 a'^ — kia^ — k^aia^ 



k' 2 a' x — k 3 a\ — kia 2 a^ 



k 2 a 3 



0,20,3 + kik 3 a 2 
= kids - fc 3 a*a 2 
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with k' = yj\ — 0202- This map defines a canonical 
transformation of the Poisson algebra, which in terms of 
the angles reads: {aifij} = {ctij, ay} = f3j} = 

0. We now canonically quantize this theory, i.e., by re- 
placing the angles by Hilbert space operators and the 
Poisson brackets by commutators, so that [a, f3] = (h, 
where £ is a complex parameter. It can be shown that 
the quantum operators corresponding to k,a,a*, satisfy 
the commutation relation of the g-oscillator algebra 

qa) a — q~ 1 aa) = q — q^ 1 



ka) = qa)k ka = q 1 ak , 



(10) 



with quantum deformation parameter q — eS h and k 2 = 
q(l — a^a). Upon this quantization, the map -R123 be- 
comes a quantum operator satisfying by construction the 
quantum tetrahedron equation ([5]). It has been shown in 
[111 ] that it is possible to construct the matrix elements 
(n[, n' 2 , n' 3 \R\ni, n 2 , n 3 ) in the basis of the Fock space 
constructed from the g-oscillator algebra. The opera- 
tor Rijk define an automorphism of the triplets of tensor 
products of the g-oscillator algebra Of 3 — > Of 3 , and it 
also has the property of being non-degenerate in F® 3 . 
This fact, together with the tetrahedron equation (O, 
implies the validity of the standard Yang-Baxter equa- 
tion (which signals the integrability in two dimensions) . 
In fact, tracing out in the Fock space F a , the following 
equations are obtained: 



Rb c RbdRcd — RcdRbdRbc 
LvbLvcRbc = RbcLvcLyb 



(11) 

(12) 



One affinization of the solution of [11J has been given in 
[T2I ] as follows: consider the layer to layer transfer ma- 
trix T mn (l,fj.) which are related to the L-operators by 

T mn (l,fJ.) = Iir=i Tlj=m L j( l i^j) and ma y be obtained 
from an ansatz that solves the local Yang-Baxter equa- 
tion: 



Li,2{u 3 ,l 3 ) 





^3 







,t 






-q~H 3 fia 3 /i 3 fc 3 
-q-H 3 fj, 3 



where we have chosen the Fock space F a as F 3 , for conve- 
nience. It has been shown that the solution of the affinc 
TE (when all parameters U are equal to each other) has 
symmetry U q (sl{n)), where n is an arbitrary integer and 
can be considered as the emergent coordinate of the third 
dimension. 

Equations (fTTj) (fP2|) can be associated either to an in- 
tegrable two-dimensional statistical system or to a quan- 
tum system in (1 + 1) dimensions. Therefore, the new 
solution of the Zamolodchikov equation tells us that is 
possible to break up the (2 + l)-dimensional lattice sys- 
tem (in a consistent fashion) into (l + l)-dimensional ones 



for each row and column of the lattice or into classical 
statistical systems in (2 + 0)-dimensions at a fixed time. 
This result relates the one-dimensional Mott point (i.e., 
the critical value of t/U in the one-dimensional analogue 
of |l| and its order parameter) to its two-dimensional 
counterpart as follows: on the one hand, in each row or 
column the XXZ spin system with hamiltonian, 



rrld 
11 XXZ 



"£(S?S* +1 + SVSV +1 AS? Sfa) + H b (13) 



where H b = a{S 7 { - S Z L ) and A = (q + q-^/2 and 
a = (q — q~ 1 )/2, is know to posses the symmetry 
U q (sl(2)) gs U q (sl(2)) Note that q = 1 at A = 1, 

so that the boundary conditions are irrelevant, and we 
can choose periodic boundary conditions without break- 
ing the quantum group symmetry. The system defined by 
(|13p can be mapped to a one-dimensional nearest neigh- 
bors Hubbard (fermionic) system (using a Jordan- Wigner 
transformation) and may also be bosonized as a Luttingcr 
system, which is a conformal field theory (CFT) [24| with 
central charge c = 1 [15| [16[ . It has been shown that this 
system undergoes a Mott transition at A = 1 be- 
coming an insulator. Morever, at the Mott transition 
(with q — 1), the quantum group symmetry becomes 
the ordinary su(2) symmetry (as far as the algebraic and 
symmetry properties are concerned, we treat su(2) and 
sl(2) as interchangeable) and the CFT will have sym- 
metry su(2) 1 eg) su(2) 1 (because two CFTs, one for each 
chirality are needed) which becomes su(2)®su(2) at long 
distances. This CFT can be realized by a Wess-Zumino- 
Witten (WZW) model at level k = 1, or by two chiral 
bosons (of opposite chirality) compactified on a circle at 
the self-dual radius (under i?-duality) , where the vertex 
operators, known as currents J* = e lv/ ^,are well defined 
under a shift in the zero-mode of the fundamental bosonic 
field of the WZW model <j> — > <fi + 2-nr (r is the compact- 
ification radius). The currents have scaling dimension 
(1, 0) and together with J 3 (z) = id(f>, of the same dimen- 
sion, satisfy the su(2) 1 algebra. The i?-duality in these 
systems is well-known ( 24]), for which an exchange of r 
for 1/r leaves the spectrum unchanged but the elemen- 
tary degrees of freedom are exchanged between charges 
and vortices. 

For A = 1 + e, the system ([TBf develops an energy 
gap in the spectrum E g = 4exp(a/e) and exhibits a 
charge density wave (CDW) order parameter [6| , defined 
by (p(i)) = 1/2[1 + where P = (VKOV'iW + 

ipl(i)ip r (i)) = 1/^/eexp (— a/y/e) (a is a constant, i la- 
bels the lattice site and the expetaction value is taken 
in the ground state). The corresponding low-energy ef- 
fective theory is given by the Sine-Gordon theory for the 
non-chiral effective bosonic field 0, of Lagrangian density: 



^sg 



-{defy 2 +J0COS-V/20 



(14) 
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The last term can be interpreted as having origin in 
the Umkplap processes naturally arising in the lat- 
tice fermionic description, in which it is given by 
(ipL^L^RipR.) (where Vl = exp(-i0 L /v / 2) and ip R = 
exp(f(ftp/\/2), with <f> — 4>l + 4> R - As it has been pointed 
out in [17], for f3 > (repulsive interactions), the ground 
state energy is minimized for a constant ground state 
density {4>) = n/y2, which in the fermionic representa- 
tion corresponds to (ip R i>L) = —{^r^l) = if, where / 
is a real function that changes sign under a tt/2 rota- 
tion around any axis (with appropriate generalizations 
for the one-dimensional case [17]]) . Therefore, this one- 
dimensional system exhibits a d-wave density order pa- 
rameter, meaning that on the ground state the quan- 
tity (ip^(k)ip(k + it /(a))) breaks individually time rever- 
sal, translation invariance by one lattice site and tt/2 
rotation (around any axis) symmetries, while preserv- 
ing the composition of any two of them. (Note that 
(ipftipL) = — (^r^l) yields P = 0, implying a constant 
density profile. 

As consequence of the projection-like character of the 
solution of [ill], by consitency, the system on the lattice 
must be in a (two-dimensional) d-density wave and also 
will be described by a CFT with c = 1 in all phases. 
In orther to further account for the charge neutrality of 
the spin excitations and the additional lattice symmetries 
H, the CFT should be taken on the orbifold S 1 /Z 2 and 
modded out by the lattice symmetry D4. This CFT has 
been identified as characterizing the critical point of the 
six- vertex model or the four-state Potts model [l8[ . 

On the other hand, the lattice fermion system can be 
viewed in an alternative way: consider a two-layered sys- 
tem periodic in the 'temporal axis', and trace out over the 
temporal dimension. The resulting Yang-Baxter equa- 
tion has symmetry U q (sl(2))®U q {sl{2)) corresponding to 
the symmetry of the six- vertex model fl3j ] . whose phase 
diagram depends on the Boltzmann weigths at each ver- 
tex of the lattice throught the parameter Agt, = (a 2 + 
b 2 - 2c 2 )/2ab, where a = exp(— /3E a ), b — exp(— /iE h ), 
and c = cxp(— f3E c ) are the weights at each vertex. The 
transfer matrix of the six- vertex model is given by that of 
the XX Z model through a Wick rotation. At the Mott 
critical point this implies that (A = — Aq v = 1) it is in 
the antiferroelectric phase. ( Note that the spins are also 
rotated changing the sign of the spin- wave term [151 ]) 

Now we would like to write down an effective field the- 
ory (EFT) for the model on the square lattice at the Mott 
critical point [3] . We first choose the effective degrees of 
freedom, and impose their characteristic symmetries on 
the theory under construction. In our case, both sym- 
metries are given by the exact solution discussed above. 
As it was first pointed out by Witten, there is a close re- 
lationship between quantum groups, vertex models and 
Chern-Simons (CS) gauge theories: the expectation val- 
ues of the Wilson loops can be calculated as statisti- 



cal sums of Boltzmann weights in suitable defined ver- 
tex models, so that the mathematical structure of quan- 
tum groups encodes the topology of planar Wilson loops 
[22! ]. However, CS theories posses naturally the sym- 
metry U q {sl{2)) (with q — exp(2ni/k) where k is the CS 

coupling constant) ^ [27j], and not U q (sl(2))®U q (sl(2)) 
. This mismatch is a consequence of the absence of par- 
ity conservation in the CS gauge theories. The simplest 
CS-type theories that preserve parity are the double CS 
gauge theories (which contain two u(l) chiral gauge fields 
of opposite chirality, namely right and left) [201 ] : 

k f k f 

Sdcs — 7- I d 3 x a R A da R - — I d 3 x a L A da L (15) 

where a R (a^) denotes the right (left) gauge field. This 
theory is known to be equivalent to the BF theory 
PH and it can also be written as a mixed CS theory. 
Here we are using a A da as a short-hand notation for 
the lattice version of the CS coupling a^K^^a^ with 



K, 



= S. 



, a ,vd a , S^f{x) = f(x + ae^), S fl f(x) = 
(f(x + ae M ) — f(x))/a, (with a the lattice spacing). At 
the Mott transition we have k = 1, since the coupling 
k fixes the unit of charge and the statistics of the exci- 
tations, and we find consistency with the fact that the 
effective degrees of freedom are density waves of the un- 
derlying electron system. Note that the u(l) CS theory 
can be considered as the broken parity phase of the su(2) 
CS theory, where the relation to the six-vertex model has 
been established [22| [23j |. 

We now impose periodic boundary conditions to the 
EFT, i.e., compactify the space domain on a torus. Cut- 
ting down the torus along any cycle, induces the loos- 
ing of gauge symmetry on the cycle, so that the gauge 
fields become boundary dynamical degrees of freedom 
[i3][3l|[33] which are free chiral bosons (c = 1 CFTs), 
representing charge density waves (CDW) described also 
by Luttinger systems . In the quantum theory obtained 
after quantizing these classical bosonic waves, there is a 
shift in the coupling parameter k that is properly taken 
into account by the Sugawara construction: k — > k + c v , 
where c v is the dual coxeter number of the symmetry 
algebra of the gauge group ( c v = 2 for su(2)). How- 
ever, the identification of the Mott transition is done at 
the classical level, implying that the topological order re- 
mains given by the relation q — exp(27ri/fc), with no shift 
in k. 

We would now like to show that in the EFT, the emer- 
gence of a d-wave order considered before is natural. For 
this, we focus on the electric-magnetic duality between 
charges and vortices implicit in the EFT. Let us consider 
the charge current j^ix) degrees of freedom in the direct 
lattice, and corresponding vortex current M (JQ) in the 
dual lattice (whose sites are in the center on each paque- 
tte of the direct lattice). We assume that these degrees of 
freedom can couple, and the low-energy action for their 
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interaction is given by a mixed Chern Simons theory 

k f k f 

Smcs = -7— I d 3 xaAda—— / d 3 x b A db + a, t /J (16) 

47T J 47T 7 

where we have introduced two gauge fields a M and 6 M 
for the current and vortex degrees of freedom, respec- 
tively. The relevant definitions are — kK^ tV a v , ^ — 
kK^ vbjji — S mu c^^ a ^i/d a7 S^f^x^) f \x -\- &€^), 

dfj,f(x) = (f(x — ae M ) — f(x))/a. In the dual lattice, 
the system dual to the original one is a two-dimensional 
XX Z spin system with coupling constant A -1 . The de- 
grees of freedom corresponding to these spin variables 
are vortices. For A < 1 the spins in the dual system are 
frozen, there are no spin waves, and the system is in the 
antiferomagnetic phase. Therefore, its effective action is 
given by a CS theory with punctures [14J: 



5, 



cs 



4-7T 



J d 3 x a" K^ax + ^ 4>° [5{x d , y d ) 



- S(x d + l,y) - S(x d ,y d + 1) + 8(x d + l,y d + 1)] , 

where a\ is a (different) abelian CS field and means 
that the sum is taken over all fundamental domains . 
Each domain has period 2a and contains four vortices 
in antiferromagnetic array. Therefore, the classical low- 
lying states reproduce an antiferromagnetic current pat- 
tern. At the quantum level, Gauss' law selects the phys- 
ical states from the lattice CS gauge theory with punc- 
tures. The quantum order of the ground state of this 
theory (staggered flux phase) breaks translation and par- 
ity symmetries by one lattice site and also time reversal 
invariance, but it is, however, invariant under the com- 
position of any two symmetry transformations, satisfy- 
ing the definition of the d-density wave order invariance 
([l3|). Therefore, the EFT analysis shows that the Mott 
transition for the system defined by (JTJ) is of topological 
(i.e., Kosterlitz-Thoulouse type) nature, characterized by 
the emergence of CS vortices in antiferromagnetic arrays. 
Note that the resulting theory breaks the Z2 vortex sym- 
metry associated and, therefore, the two-dimensional chi- 
rality on each plaquette (fTS)) . 

Finally, we would like to discuss the behavior of the 
previously considered EFT away from the Mott critical 
point. By analogy with the CS theory of the quantum 
Hall effect, we could expect a ground state stable against 
small doping. In that case, for the simplest Laughlin in- 
verse filling fractions k = m (m odd integer), the ground 
state is a droplet of incompressible quantum liquid [32| 
(however, other phases with more exotic quantum orders, 
like Nematic phases are also possible in other regimes) 
and is stable under small perturbations away from the 
center of a given plateu in the conductivity. In the Mott 
system, we have already assumed that the dynamically 
generated vortices act as external statistical fields for the 
new electrons injected in the system by doping (this can 



be considered as an extension of the i?-duality).At the 
self-dual point, statistical magnetic fields can be inter- 
changed with statistical electric fields (on a torus). After 
imposing the lattice symmetries, the low-lying effective 
Hamiltonian for the injected electrons (in first quantiza- 
tion) is: 



^-E^ 2 i(^ + ^) + ^- yi 2 )], (17) 



where can take the values ±A Therefore, the electric 
potential changes sign in x — ±y, producing domain 
walls between regions with different electron densities. 
Similar results can be obtained using the W4 symme- 
try, which is related to the relevant perturbations of the 
Ashkin- Teller and six-vertex models away from the crit- 
ical point [m . 

One consequence of having discussed the EFT is that, 
a-posteriori, the behavior of the electrons can be more 
easily understood. It can be shown that the interaction 
term in the hamiltonian {T]) in the the continuum limit 
contains a chemical potential term of the form —fx p, 
with fj, = A , which ensures the half-filling condition. 
Therefore, changing the chemical potencial by doping in 
Sfj, modifies the hamiltonian (in the spin representation 
© by H(A) H(A) + Sfij: {ij) S?S]. At A = l,the 
dynamics of the electron system is given by the dou- 
ble CS theory (fT5j). whose hamiltonian can be defined 
as the temporal component of the stress-energy tensor 
H cs = Too, where = SS^/Sg^ and where is the 
metric tensor. However, the CS action is topological and, 
therefore, independent of the metric and H cs = for each 
chiral component leading to H(A = 1) = 0. This means 
that doping the system away from the critical point, the 
dynamics is controlled by an effective Ising hamiltonian. 

To sum up, our study of the EFT for the model ([I]) 
shows that the Mott transition is of topological (i.e, 
Kosterlitz-Thoulouse type) nature, characterized by the 
emergence of CS vortices in antiferromagnetic arrays. 
The symmetry su(2) 1 ®su(2) 1 of the critical point implies 
that electric and magnetic vortices can be intercanged, 
providing effective atractive and repulsive forces. Doping 
the system produces domain walls, signaling a quantum 
gas-liquid phase coexistence, which belongs to the Ising 
universality class. 
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